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We present a new method of obtaining nonlinear integral equations characterizing the thermo- 
dynamics of one-dimensional multi-component gases interacting via a delta-function potential. In 
the case of the repulsive two-component Bose gas we obtain a simple system of two non-linear in- 
tegral equations, allowing for an efficient numerical implementation, in contrast with the infinite 
number of coupled equations obtained by employing the Thermodynamic Bethe Ansatz. Our tech- 
nique makes use of the Quantum Transfer Matrix and the fact that, in a certain continuum limit, 
multi-component gases can be obtained from appropriate anisotropic spin chains. 
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Introduction. Recent advances in magnetic and opti- 
cal trapping of ultracold quantum gases have opened new 
possibilities in investigating the physics of strongly inter- 
acting particles in one-dimension Coupled vifith the 
fact that the strength of the atomic interaction can be 
controlled using magnetic Fesbach resonances or state- 
dependent dressed potentials Jj, these new experiments 
provide a fertile ground for testing the theoretical predic- 
tions obtained from the study of various integrable mod- 
els solvable by Bethe Ansatz. As a result, we have wit- 
nessed the experimental realization of the Lieb-Liniger 
model [3" using bosonic atoms [IHS] and, moreover, 
the measured thermodynamics (Hj is described very well 
by the Yang- Yang thermodynamics, which was developed 
in [lOj . using, what we call nowadays, the Thermody- 
namic Bethe Ansatz (TB A) [11] . In the case of integrable 
multi-component ID gases, TBA produces an infinite 
number of integral equations [TTl [T2J , which makes the 
extraction of physical information and comparison with 
prospective experimental data difficult. In this Rapid 
Communication, we propose a solution to this problem 
using the Quantum Transfer Matrix [T31II1], and the con- 
nection between multi-component gases and anisotropic 
spin-chains. Our method has the advantage of providing 
a finite number of non- linear integral equations (NLIE), 
which are extremely suited for numerical computations, 
in stark contrast with the TBA result. 

Model. We consider a one-dimensional system of M 
bosons, of equal mass m = 1/2, with n internal degrees of 
freedom, interacting via a delta-function potential. The 
many-body Hamiltonian is 

i—1 * i<j i — 1 

where c is the coupling constant and we consider h — 
kg — I. The first and the second term in the Hamilto- 
nian ([T]), represent the kinetic and and the interaction 
energy, while the third is the Zeeman term, where Mi 
is the number of particles in the hyperfine state \i) and 
fii the respective chemical potential. The interaction is 



attractive for c < and repulsive for c > 0. The model 
is integrable and was solved in [31 [TSHTT] . 

Despite the fact that the model is integrable, com- 
puting the thermodynamics is still an incredibly difficult 
task. The first attempt was done by Yang and Yang 
[10) for the spinless bosons using the Thermodynamic 
Bethe Ansatz. TBA has proved to be an extremely useful 
technique and was used to obtain information about the 
thermodynamics of various exactly solvable models [TT] , 
however, it has one big shortcoming: for a large class of 
models it produces an infinite number of coupled NLIE. 
This is also the case for our model, with the exception of 
the spinless case {n = 1). From a very simplified point 
of view this is due to the fact that the spectrum of the 
theory contains infinitely many branches. The numeri- 
cal implementation of this system of equations requires 
various truncations, which introduce an uncontrollable 
source of numerical errors and makes the extraction of 
relevant physical information extremely difficult. There- 
fore, it is highly desirable to devise a numerically efficient 
method, which provides a finite number of NLIE. 

Fortunately, in the case of lattice models such a 
method exists. Developed in [131 dH) the Quantum 
Transfer Matrix (QTM) approach was successfully ap- 
plied to a large class of integrable spin-chains and even 
models of strongly correlated electrons. Within this ap- 
proach, the thermodynamic properties of the model are 
obtained from the largest eigenvalue of the QTM, which 
in the thermodynamic limit gives the free energy. Even 
though, there is no equivalent of the QTM for contin- 
uum models, it is well known [THl [H] , that a large class 
of integrable models can be obtained from lattice ones 
in suitable continuum limits. We say that an integrable 
model is the continuum limit of a lattice integrable model 
if, by performing this limit in the Bethe equations and en- 
ergy spectrum of the lattice model one obtains the Bethe 
equations and energy spectrum of the continuum model. 
The natural consequence is that the thermodynamics of 
the continuum model can be obtained from the thermo- 
dynamics of the lattice model, if we take the same limit. 
This method will be used in this Rapid Communication, 
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allowing us to obtain an efficient thermodynamic descrip- 
tion of the two-component Bose gas (for the single com- 
ponent case see [10]). 

In this publication we study and give results for the 
Hamiltonian ([T]) in dimensionless units, however physical 
units can be restored easily. For particles with mass m 
and contact interaction strength 5, see e.g. (1) in |1T], the 
units of temperature, chemical potential, magnetic field, 
particle density and susceptibility, compressibility, heat 
capacity, and entropy per length are Tq = / {2ma?kB), 
Mo = ^0 = fi^/(2TOa^), do ~ 1/a, Xo — 1^0 = 2ma/h?, 
cq — So = ks/a- The quantity a is a length scale that 
can be chosen freely yielding the dimensionless coupling 
constant c = mga/h^ appearing in (1). In all figures 
presented in this paper, physical data are shown in the 
given units and for dimensionless coupling c = 1 which 
is realized for any parameter values of m and g with a 
suitably chosen a = ti^/{mg). 

Method. The general strategy is the following: first, 
we identify the spin-chain from which we will obtain in 
the appropriate continuum limit our model of interest. 
The next step is the investigation of the thermodynamics 
of the lattice model using the QTM technique, which 
will provide a finite number of NLIE. Finally, the desired 
result is obtained by taking the continuum limit in the 
equations for the lattice model. Let us make these general 
considerations more precise. In the case of the repulsive 
two-component Bose gas, the relevant lattice model is the 
Uq{sl{3)) Perk-Schultz spin-chain P^VES^, with periodic 
boundary conditions characterized by the Hamiltonian 



(2) 
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a, 6=1 



isin7 sign(a-6)ei^eif ^')~5]^/i,ei^. 



L 3 



a.b=l 



j = l a=l 



where 7 G (0,7r) determines the anisotropy, (g = e*'*'), 
L is the number of lattice sites, J > fixes the en- 
ergy scale, ha are chemical potentials and e^^jj = (g) 
Gob ® /"^^"J , with Cab the 3 by 3 matrix with elements 
{fiab)ij — Sai5i,j. This Hamiltonian is the sum of two 
terms. Hi, (in the brackets), and the chemical poten- 
tial term He, which does not break the integrability 
of the model (24]. The Uq{sl{3)) Perk-Schuhz spin- 
chain can be obtained as the fundamental spin model, 
(see ^7\), associated with the trigonometric Perk-Schultz 

R-matrix |M|, ^iv) = El.b,c,d=iKdi^Kb e,d, with 
Kliv) = '-^^,Kbiv) = 'i^, (« ^ b),Rtliv) = 
e^signia-b)v ^ (^a ^ b), and all the other elements zero. 
The transfer matrix T(t;) is the 3^ by 3^ matrix with el- 
ements Tb(w) = E{c} Uj=i K]*'b-iv) ' where a,b, and 
c are multiple indices, i.e., a — (ai, • • • , a^) and the sum 
is over all c with ci = c^+i. The Hi, part of the Hamil- 
tonian ([2]) can be obtained as the logarithmic derivative 



of the transfer matrix Hj, = J sin 7 \ y^o. In order 

to treat the model at finite temperature we introduce the 
3^ by 3^ Quantum Transfer Matrix with elements 

N/2 

{tqtm)Uv) = E ''''' n Rcr"4:- + 

{c} J = l 

where Ci = Cn+i, Kdi'") = Rda(«): " = -Jpsinj/N, 
with N the Trotter number and f3 ~ 1/T. The largest 
eigenvalue of the QTM, from which the free energy of 
the model can be obtained, has the form: Aqta/(i^) — 
E-=iA,(«) I2gj,with 

X / \ ^ ^ \j. ^ .<lj-iiv-ij)qj{v + i-f) gf^. 



where 



0±(«) = ( 



sinh(D±m) 



N/2 



and 



nfi^inM'^-^r^') J = 1,2 
(l)+{v) i = 3. 

The two sets of parameters, {wr"*^}, are the Bethe 
roots of the QTM and satisfy the Bethe equations 
Xj{vi^^)/Xj+i{vli'^) = -1, r = I,-- - ,7V/2. The Bethe 
roots are distributed in certain strips of the complex 
plane which are independent of temperature and Trot- 
ter number. This type of specific distribution allows 
for the definition of auxiliary functions and the use of 
Cauchy's theorem in deriving NLIE for the largest eigen- 
value [m results, which in the contin- 
uum limit produce Eqs. [3] and [4j can be understood as 
the mult i- component generalization of the equations pre- 
sented in [22]. 

Main result. In order to obtain the thermodynam- 
ics of the two-component Bose gas we have to perform 
a particular continuum limit in the NLIE obtained for 
the Uq{sl{3)) spin-chain. The spin-chain is characterized 
by the following parameters: lattice constant S, num- 
ber of lattice sites L, anisotropy 7, strength of the in- 
teraction J, and chemical potentials /ii, /12, ^3- The two- 
component Bose gas is obtained by performing the limit 
7 = TT - e , (5 ^ 0(e2) , L 0(l/e2) , J 0(l/e4) , hi 
0(l/e^) with e — >■ 0. Performing this limit in the Bethe 
equations and energy spectrum of the Uq{sl{3)) spin- 
chain, we obtain the Bethe equations and energy spec- 
trum of a two-component Bose gas characterized by the 
parameters: length I — LS, coupling constant c = e^/(5, 
mass of the particles m — JS^, and chemical potentials 
/ii = Je^ + h2 — hi,^2 — h3~h2. In order to simplify the 
formulae, we are going to consider JS^ = 1 (mass of the 
particles equal to 1/2) and introduce ^ = (/ii+/i2)/2 and 
an effective magnetic field h = (/ii — /i2)/2, where fix 2 
are the chemical potentials of the spin up and down par- 
ticles. We are now ready to state the main result of this 
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Rapid Communication. The thermodynamics of the two- 
component repulsive Bose gas is completely characterized 
by the following system of nonlinear integral equations: 

Inai(fc) = -I3{k^ - ^-h) 

+ [Kq * In Ai]{k) + [Ki * In A2] {k + ^e) , (3a) 
lna2(fc) = -/3(fc2 + 

+ [K2 *\nAi]{k- le) + [Ko * In A2] {k) , (3b) 

where A,{k) ^ l + a,{k) , Ko{k) = 2c/{P + c^) , Ki{k) = 
c/[k{k + ic)],K2ik) ^ c/[k{k - ic)] and [/ * g\{k) = 

^ J-^ f{k — k')g{k') dk'. The grandcanonical potential 
per length is given by 



T 
'2t: 



-\-oo 



{In Ai{k) + \n A2{k)) dk, (4) 



from which all the other thermodynamic quantities can 
be obtained. Eqs.[3]and|4]are new in the literature, and 
to our knowledge this is the first example of an efficient 
thermodynamic description for a multi-component con- 
tinuum integrable model at all values of the relevant pa- 
rameters (temperature, chemical potentials and coupling 
constant). 




FIG. 1. (Color online) Population densities di, d2 in tlie up- 
per, lower part of the figure as functions of temperature T for 
jj, = 15, c = 1 and various effective magnetic fields h (in units 
of do, To, Ho and ho). 

Analytic limits. First, we will consider the limit 
c — > 0. Then, we have limc-^o c/(A:^ + c^) = T:5{k) and 
limc-j-o c/[fc(fc ± ic)] = and the integral equations de- 
couple 

Inai(fc) = -I3{k^ ^ h - h) + ln(l + ai{k)) , 
lno2(fc) = -I3{k^ -fi + h)+ ln(l + a2{k)) . 

These equations are easily solvable, and we find 

rp /'+00 




FIG. 2. (Color online) Susceptibility x as a function of tem- 
perature T for fi — 15, c = 1 and various effective magnetic 
fields h (in units of xo, To, Mo and ho). Inset: polarization 
as a function of temperature for the same parameters. In the 
absence of the magnetic field the polarization is zero. 



which is exactly the grandcanonical potential of two non- 
interacting bosonic gases at different chemical potentials. 

In the case of an extremely strong magnetic field, h — > 
00, we expect that the system will become fully polarized 
and we should obtain the Yang- Yang thermodynamics 
[TIJ| of the repulsive one-component Bose gas. In this 
limit a2{k) ^ and we obtain 



Inai(fc) = -P{k^ 



h) + [Ko*Hl + ai)]{k)^ 



and 



27r 



In 



(1 - e^'^C' )(1 - e-'9('^ ^^+'')) 



dk . 



— ^ /_oo ^^^^ ~^ ai{k))dk , which is the result 
obtained in [10]. The same result is obtained in the low 
temperature limit (T <C /i, h, c), when the magnetic field 
is finite and fixed which shows the ferromagnetic (fully 
polarized) nature of the groundstate. 

As we have mentioned, the numerical implementation 
of the infinite number of coupled NLIE obtained with 
the help of TBA [H] is extremely difficult and encoun- 
ters serious problems in the regime of low magnetic field 
(h <C T, /i, c), or low temperature limit [Q < T <^ h, ^,c). 
In contrast, Eqs. [3] are easily implementable and pro- 
vide reliable results in a larger domain of the parameters 
space. We have checked our results with available nu- 
merical data obtained by using the much more involved 
TBA equations J33] and found perfect agreement. In our 
case, the regime in which the accuracy is decreasing, is 
given by the low temperature and low magnetic field limit 
(T — ^ 0,h 0). This is a consequence of a first or- 
der phase transition resp. coexistence of fully polarized 
phases at T and h equal to 0. It can be seen clearly in 
Figs. [1] and [2] that, in the absence of the magnetic field the 
population levels, = - + (-l)»-i||) /2, z = 1,2, 

are equal and, consequently, the polarization defined as 
P = [di — d2)/d with d — di-\' d2 is zero. In the presence 
of the magnetic field the groundstate is ferromagnetic. 
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FIG. 3. (Color online) Specific heat C as a function of tem- 
perature T for /i = 15, c = 1 and various effective magnetic 
fields h. Inset: entropy S = — |^ as a function of tempera- 
ture for the same parameters. (AH quantities in units of co, 
To, /^o, ho and So.) 



therefore, the polarization at T = is -|-1, for any posi- 
tive field (—1 for negative field) and it decreases at higher 
temperatures. The ferromagnetic nature of the ground- 
state, which can also be seen in the T^/^ behavior of the 
specific heat C — at zero magnetic field (see Fig.|3]), 
is in accordance with a more general theorem of Eisen- 



berg and Lieb [31] on systems with spin-independent in- 
teractions (see also [2l])- Another consequence of the 
phase transition is the divergence of the zero-field sus- 
ceptibility, X = ^ feature which can be seen in 
Fig. [2] The specific heat, magnetic susceptibility, and 



compresibility, k 



(see Fig. 



monotonic behavior with local maxima 



4), present a non- 



shifting to higher 



temperatures as the magnetic field increases. 

In conclusion, we have presented a new method of 
obtaining a finite number of NLIE characterizing the 
thermodynamics of integrable multi-component ID Bose 
gases, which allows for an efficient numerical treatment 
and has significant advantages over the TBA result. 
Fermionic gases can be treated in a similar fashion, but in 
this case the relevant lattice model for the n-component 
system is the Uq{sl{n\l)) Perk-Schultz spin-chain. The 
derivation and a detailed analysis of Eqs. [3]and|4]will be 
presented in a future publication. 
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FIG. 4. (Color online) Compressibility function of tem- 
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